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Abstract. We have investigated the 2-D stellar distribution in the outer parts of three nearby open clusters: NGC 2287 
(= M41), NGC 25 16, and NGC 2548 (= M48). Wide-field star counts have been performed in two colours on pairs of digitized 
ESO and SRC Schmidt plates, allowing us to select likely cluster members in the colour-magnitude diagrams. Cluster tidal 
extensions were emphasized using a wavelet transform. Taking into account observational biases, namely the galaxy clustering 
and differential extinction in the Galaxy, we have associated these stellar overdensities with real open cluster structures stretched 
by the galactic gravitational field. As predicted by theory and simulations, and despite observational limitations, we detected a 
general elongated (prolate) shape in a direction parallel to the galactic Plane, combined with tidal tails extended perpendicularly 
to it. This geometry is due both to the static galactic tidal field and the heating up of the stellar system when crossing the Disk. 
The time varying tidal field will deeply affect the cluster dynamical evolution, and we emphasize the importance of adiabatic 
heating during the Disk-shocking. In the case of NGC 2548, our dating of the last shocking with the Plane (based on a tidal 
clump) is consistent with its velocity. During the 10-20 Z-oscillations experienced by a cluster before its dissolution in the 
Galaxy, crossings through the galactic Disk contribute to at least 15% of the total mass loss. Using recent age estimations 
published for open clusters, we find a destruction time-scale of about 600 Myr for clusters in the solar neighbourhood. 

Key words, open clusters and associations: general - open clusters and associations: individual: NGC 2287, NGC 2516, 
NGC 2548 - stars: Hertzsprung-Russel (HR) and C-M diagrams - Galaxy: kinematics and dynamics 



£h ! 1- Introduction 



Galactic open clusters (OCs) have long been shown to be ex- 
tremely valuable laboratories for many domains of astronomy. 
In particular, observations of their spatial structures can help to 
constrain dynamical models and A^-body simulations of small 
(N rs 10 3 ) stellar systems. Extended over some parsecs, OCs 
are made up from many dozens to a few thousands stars. Their 
sparseness makes them short-lived stellar concentrations. 

In this manner, most of galactic OCs evaporate entirely in 
some 10 8 years: on the 1200 or so objects of this type in the 
catalogue of open cluster data compiled by Lynga (1987), only 
about 70 are known to be older than 1 Gyr. Indeed, as a result 
of 2-body relaxation, several stars can acquire positive energy 
during the strongest (nearest) interactions with other members, 
and then leave the cluster which slowly vanishes. Stellar evo- 
lution (mass loss) and encounters with giant molecular clouds 



(GMCs, see Wielen 1991 ) also contribute to noticeably reduce 



the stellar system lifetime. In addition, the galactic potential 
plays a major role in the disruption process by its regular grav- 
itational harassment, through tidal disruption. 

As a consequence, only the initially richest clusters (which 
are more gravitationally bound) and those situated at large 
galactic radii (where the probability of encounter wi th a G MC 
is lower) can live a few Gyr. De la Fuente Marcos ( 1998 ) has 
shown that our Galaxy should host several hundreds of thou- 
sands of open cluster remnants which have been disrupted by 
the Galaxy and by the ir own dynamical evolution (evapora- 
tion). Bica et al. ( 2001 ) indeed observe several dissolving star 
cluster candidates. 

Orbits of OCs being quasi-circular, with Z-oscillations of 
small (< 1 kpc) amplitude, their trajectory includes many pas- 
sages through the Disk: each of these gravitational shocks heat 
up and compress the cluster which then takes a prolate shape 
flattened to maximum at Z = (Leon 1998 1. Repeated disk- 



Send offprint requests to: G. Bergond (gbergondOeso . org) 

* Plate scanning was done at the Centre d' Analyse des Images 
(CAI) with M.AM. A. (Machine Automatique a Mesurer pour 
l'Astronomie), a facility located at the Observatoire de Paris, de- 
veloped and operated by I NSU (Institut National des Sciences de 
1'Univers, CNRS). Web site http://dsmama.obspm.fr 



shockings speed up the di sruption, as studied in globular clus- 
ters (Combes et al. 1999) and after each passage through the 



Plane, members rejected in the halo of the system are stripped 
out by the gravitational field of the whole Galaxy. These ejected 
stars will form a tidal tail which extends far from the cluster in- 
ner regions. Structural studies of the surroundings of OCs can 
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thus provide much information on the past and present dynam- 
ical evolution of the cluster. 



Terlevich (1987) and de la Fuente Marcos (1997 and ref- 
erences therein) performed realistic simulations of open clus- 
ter evolution. From the observational point of view, Grillmair 



et al. (1995) and Leon et al. (200C) detected important tidal 



tails extensions around galactic globular clusters using a wide- 



field star count analysis. Leon et al. (1999) also found notice- 
able tidal tails around several binary star clusters in the Large 
Magellanic Cloud. 

The main aim of this study is to detect such tidal tail sig- 
natures in the overall cluster structure around selected galac- 
tic OCs, by performing multicolour, wide-field star counts 
on Schmidt pho tograp hic plates. First results are presented in 
Bergond et al. (2001). In Sect. 2 we discuss the selection of 
three open clusters. Section 3 presents the data analysis, and 
a brief discussion on our observational limitations. In Sect. 4 
we focus on each of the three observed OCs. Next, in Sect. 5, 
we discuss these results together with numerical simulations, 
and their general implications for galactic open clusters, before 
summarizing this work in the last section. 



2. Open cluster sample 

It is expected that tidal tails are very weak stellar overdensi- 
ties, and to better detect them it is necessary to reject at max- 
imum areas with strong background variations due to a differ- 
ential interstellar extinction. A hole of extinction can indeed 
artificially enhance the star-count in this more transparent area, 
giving the impression of a higher density of stars which we 
could mistake for a tail. For this reason, we have made a first 
selection, restricting our choi ce to the few OCs which present 
b > 10° (see Janes & Adler 1982 for the galactic distribution 
of OCs). Such high latitudes allow us to avoid the major part 
of dust clouds which are generally well concentrated along the 
Plane Z = 0, w ith a mean scale height of 160±20 pc (Pandey 
&Mahra 119871). 



However, the extinction within the Plane is highly variable 
from one direction to another, and there are windows of trans- 
parency, especially between £ = 210° and £ = 240° (Chen et 
al.~ 



1998 (. The whole map of dust emission (as traced by IRAS 
100 i-im) in these directions is visible in Fig. [l| wh ere we have 



overplotted all the known OCs in the Lynga (1987) catalogue. 

We had yet to select relatively rich clusters in order to per- 
form a reliable star-count analysis. Very poor OCs do not offer 
a sufficient contrast with respect to the field stars, and may be 
only asterisms, that is apparent concentrations of stars due to 
a pure perspective effect (e.g., Baumgardt 1998| ). At last, the 
cluster must have the appropriate dimensions: too small clus- 
ters are difficult to treat with the typical scale on Schmidt plates 
(67"l/mm) due to limitations in resolution of emulsions and 
diffusion problems in crowded areas, whereas really too wide 



clusters like Melotte 111 (Odenkirchen et al. 1998) cannot be 
studied on one single w 30 deg 2 plate. 

We finally selected as the most interesting candidates 3 
rich, bright and low-absorbed OCs: NGC 2287 (= M41), 
NGC 2516, and NGC 2548 (= M48). Table presents the pairs 



270 260 250 240 230 

Galactic longitude idcg) 

Fig. 1. Skyview map of IRAS dust emission at 100 /vm 
in the regions going from (£,b) = (222°,-28°) to (£,b) = 
(278°,+28°). All the known open clusters in the Lynga cata- 
logue are shown as small crosses (+). The three selected objects 
are indicated. They avoid strong and differential extinction ar- 
eas. 

Table 1. Scanned Schmidt plates encompassing the clusters. 



Cluster 


Plate 


Emuls. 


Filter 


Exp. 


Epoch 


NGC 


ID# 


Kodak 


used 


min 




2287 


SRC557 


IllaJ 


GG395 


60 


1979.898 


2287 


ES0557 


IIIaF 


RG630 


60 


1985.930 


2516 


E11031 


IllaJ 


GG385 


15 


1994.117 


2516 


El 1190 


IIIaF 


RG630 


15 


1994.444 


2548 


SRC775 


IllaJ 


GG395 


60 


1986.246 


2548 


SRC775 


IIIaF 


OG590 


70 


1985.057 


2548 


SRC776 


IllaJ 


GG395 


60 


1983.192 


2548 


SRC776 


IIIaF 


OG590 


70 


1985.953 



of photographic plates we scanned for this study and Table g 
summarizes the main properties of these OCs. 

3. Data analysis 

All the plates were digitized with a resolution of 10 /jm with the 
Machine Automatique a Mesurer pour VAstronomie (MAMA, 
see Berger et al. 1991). Source extraction was done using the 
S Extractor software (Bertin & Arnouts 1996), at a 3a-level for 



each pair of plates. Once the astrometry was performed in each 
colour, local software was used for the cross-correlation of the 
B and R catalogues. 

This allows us to build a colour-magnitude diagram 
(CMD), in instrumental magnitudes, for the whole plate . Next , 
we adopt the method described by Grillmair et al. ( 1995 ), 
which consists of selecting regions in the CMD where the 
"cluster vs. field contrast" is important. For this purpose, we 
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Table 2. Main characteristics of the three studied open clusters as extracted from the Base des Amas (see Mermilliod 1995) 
http : //obswww . unige . ch/webda/ , except for the masses (particularly uncertain), which come from Bruch & Sanders (1983 ) 
for NGC 2287, or from Pandey et al. (1987) for NGC 2516. Some values showing important discrepancies in the literature are 
specified in the note below the table. 



Cluster 


a 


5 


£ 


b 


Dist 


Z 


Rad. 


Mass 


Age 


Spect.-type 


[Fe/H] 


Eb-v 


name 


J2000.0 


J2000.0 


ii 


II 


pc 


pc 


pc 


M 


Myr 


on turn-off 


dex* 


mean 


NGC 2287 


06?46! n 9 


-20° 44' 


231.° 10 


-10P23 


693 


-123 


4.1 


294 


243 


B5 


+0.04 


0.03 


NGC 2516 


07! 1 58f n 3 


-60°52' 


273 ? 94 


-15P88 


409 


-112 


1.9 


170 


113+ 


B3 


+0.06 


0.10 


NGC 2548 


os! 1 ^! 1 ^ 


-05°48' 


227 ? 92 


+15P37 


769 


+204 


4.8 


N/A 


360 


AO 


+0.08 


0.03 



T Meynet et al. (|1993j) fit NGC 2516 w ith an older isochrone at 140 Myr. 
* For NGC 2516, Jeffries et al. (|l997|) proposed a value of [Fe/H] ~ -0.32. 



use a so-called colour-magnitude S/N function s(i,j) defined 
by subdividing the CMD in multiple cells (50 x 50 bins here, 
i.e. about 0.1 mag both in B and B —R). We refer to the papers 
of Leon et al. (1999, 2000| ) for details of the application of the 
method. We then retain only cells where the number of cluster 
members with respect to the field stars is the highest, by defin- 
ing a threshold we choose as a compromise between significant 
star-count (in order to prevent background fluctuations) and re- 
duced contamination by field stars. 

In the case of relatively nearby clusters, it is possible to use 
a direct "cut" in magnitude in the whole star count, selecting 
only the brightest sources. Most of weak field stars are then 
directly eliminated, whereas only the bottom of the OC main 
sequence (that is, most of late-type dwarf members) is missing. 
Keeping only bright OC members to minimize pollution also 
obviously decreases the number of detected stripped stars from 
the cluster, which are preferential ly low mass stars, as noted 
both from observations (Leon et al. 2000) and simulations (e.g., 
Combes et al. h999|), but the S/N of the tidal tails relative to the 



background/foreground stars will be higher due to the richness 
of the three selected clusters. 

The three OCs we have chosen being nearby systems (with 
a distance to the Sun in the range J— § kpc), we have applied 
cuts in instrumental magnitudes corresponding on average to 
about My ~ 5. Using the "universal" mass function of Kroupa 
( 2001 ), this implies that we may miss about 90% of the theo- 
retical cluster total mass (see Fig. [|). Note however that Scalo 
( 1998 ) claimed that a universal mass function is not justified 
empirically for stars with masses less than 1 Jfo. Recent the- 
oretical studies (Lastennet & Valls-Gabaud 1999) also tend to 



show that open cluster IMFs present some variations (which 
may be due to mass segregation effects, see the discussion in 
Sect. 5) below 1.1-1.2 Mq. 

Therefore, estimates of the fraction of the cluster mass 
we observed are accordingly subject to non-negligible varia- 
tions, and the Kroupa (2001) IMF may rather be an average of 
possibly locally different IMFs, as already observed even for 
intermediate-mass stars in OCs (e.g., Phelps & Janes 1993 or 
Sanner & Geffert |2001l). 



After this first selection from the mask in the CMD diagram 
(for example, see Fig. |J) which on average keeps N > 1500 
likely members on the whole field, we need to pay attention to 
a remaining background (Milky Way) gradient in the wide-field 



Kroupa (2001) IMF 



M (solar musses) 



Fig. 2. Kroupa ( [200 1| ) galactic initial mass function with the 
low-mass cuts (~ 1 ±0.1?Wq) applied for the three clusters. 
Clearly, we miss the numerous weak stars which populate pre- 
dominantly the tidal tails, observing possibly only 10% of the 
total cluster mass. 



encompassing the OC. Following Grillmair et al. (1995), as a 



first step we mask the cluster, up to 1 or 2 times its estimated 
tidal radius r t . The hidden cluster is then replaced by the mean 
density value just outside it (i.e., from 1.5 to 2.5r t ). The back- 
ground z is next smoothed using a fit by a bivariate polynomial 
surface: 



z(x,y) 



•ijxy 



<ij< 2. 



(1) 



A simple lxl bivariate polynomial is generally sufficient, 
except near the Galactic Plane - as it is often the case for OCs - 
where a degree 2 is preferable in order to follow the exponential 
increase of field stars towards galactic latitude zero. A higher 
order polynomial is not applicable, as it may erase similar scale 
variations such as the tidal tails we are seeking. Then, we obtain 
a raw map of overdensities 77afep> aw by subtracting the Milky 
Way gradient z from the count of likely members of the cluster: 



Tide Rav ,(x 1 y) = cmd(x,y) -z{x,y), 



(2) 
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Table 3. Parallaxes and proper motions for the three selected 
OCs, as extracted from Baumgardt et al. ( 2000| ). AT* is the 
number of HlPPARCOS members used to compute K and ^ a * s- 
fj a * includes the multiplication by the cosine of the declination. 



Cluster 


% (mas) 




^a* (mas/yr) 


^8 (mas/yr) 


NGC2287 


1.93±0.49 


8 


-4.34±0.40 


-0.09±0.40 


NGC2516 


2.77±0.25 


13 


-4.08±0.27 


-10.98±0.24 


NGC2548 


1.63±0.79 


5 


-0.50±0.70 


+0.93±0.65 



where cmd is the surface density of CMD-selected stars. 

A t this point, we use the wavelet analysis (see Leon et al. 
199^ , |2000| ). This method makes use of the so-called "a trous" 



algorithm (Bijaoui 1991 ). Relatively fast, it decomposes by 
a discrete wavelet transform (WT) the raw image in several 
planes Wj(x,y) where i is the scale of the plane, using a B-spline 
function. This kernel function Bs(x,y) is used for the recursive 
convolutions of the initial image. If we note as co the raw image 
of overdensities Tiden avl just obtained, we have successively: 



c (x,y) = Tide Rliv ,(x,y), 
x y 

Ci[x,y) = q_i xBs^,-.); 
Wi(x,y) = a(x,y)-Ci-i(x,y). 



(3) 
(4) 
(5) 



Each of the planes W, (x,y) corresponds to a certain image res- 
olution: while orders like 5 or 6 retain only "out of focus" 
features, the lower orders give the finest details. In fact, the 
last plane W; milx (x,y) - often known as LSP for Last Smoothed 
Plane - is not in the literal sense a true wavelet plane, but con- 
tains the residuals of the last convolution. However, we will 
consider it like the other planes, giving a total number of 7 
scales. This value has proved to be a good compromise between 
computing time and the maximal resolution we need. 

Finally, having obtained all the wavelet planes Wi(x,y), we 
can construct the best map to recover overdensities around the 
cluster by adding several wavelet planes. We have chosen to 
select only large scale planes related to the tidal tails sizes: the 
best results are obtained using planes from order 3 to the LSP 
(i.e., order 7 in our decomposition). The final map of tidal tails 
will then be defined as in Leon et al. ( 2000| ) by the sum: 



Tide(x,y) = £W)(x,y). 

i=3 



(6) 



This choice was also made in order to obtain a good compro- 
mise between resolution and background noise. 

4. Results 

4.1. NGC2287 (= M41) 

This bright, intermediate-age OC was particularly studied for 
prope r motions (Ianna et al. 1987 ), radial velocities (Amieux 



1988J) and spectrophotometric properties (Harris et al. |1993[ ). 
For the present wide-field structural study of the cluster, we 
found that a cut at an instrumental (R) mag. of 14 offered 
the maximum cluster contrast with respect to the field. Using 



a plate calibration based on the GSPC-II, it is equivalent to 
Bern ^ 16.2. Given a distance to NGC 2287 of about 700 pc 
and the low absorption in this direction (see Table ^|), the cut 
corresponds to main sequence stars of about My ~ 6.1, hence 
masses of the order of 0.9 Mq. The main sequence brightest 
selected stars have V ~ 9, i.e. mass greater than 3 5Wq- Taking 



the IMF of Kroupa (2001 ) within these limits, our selection of 
members may represent only « 10% of the total cluster mass. 

Despite this limitation, the overall structure observed for 
NGC 2287 (cf. Fig. |], left) corresponds relatively well to what 
is expected: a global elongated shape nearly parallel to the 
Galactic Plane (i.e., from SSE to NNW) with projected axis 
ratios estimated to 1.5:1 .0, and with extensions perpendicular 
to it (the arrow points perpendicularly towards Z = 0, a clump 
of density is visible in this direction). The dashed arrow indi- 
cates the proper motion (PM, see Tab. ||) of the OC, and the 
dashed circle traces the cluster extension (radius estimated to 
~ 35') as given in Table |[ 

The "pollution map" (see Fig. [|, right) shows no obvious 
clusters of galaxies and no anticorrelation between the tidal 
tails and the IRAS 100 map (dust absorption) which could 
bias the source count. Otherwise, it is clearly visible in Fig. |] 
(left) that the star count peak is shifted by about 15' SW from 
the center of the field, i.e. the "classical" coordinates of NGC 
2287 given in catalogues are incorrect. 

4.2. NGC 2516 



Eggen (1974) detected several members up to ^5 pc from 
the cluster center and described NGC 2516 as the Southern 
Pleiades. Indeed, this very bright and rich cluster shares a com- 
mon age and V absolute galactic motion with M45. NGC 2516 
contains many variable stars and is ideal to look for low-mass 
and active stars. It is also an essential cluster in the calibration 
of the IFMR (Initial Final Mass Relation) as it contains several 
white dwarfs. 

Among the numerous studies of this spectacular OC, the 
one by Dachs & Kabus ( |1989| ) is complete up to M v = 5!"5, 
over a relatively wide-field. They propose a total cluster mass 
of about w 1OOO5W0. 

An efficient cut at an instrumental magnitude of 13.5 
was done, as th e clus ter is particularly nearby (~350 pc, see 
Robichon et al. 1999). After cross-identification with the cat- 



alog of Dachs & Kabus ( 1989 ), this cut corresponds to main 
sequence stars brighter than V = 12.4 or My ~ 4.4. In terms of 
mass, the cut retains only stars with M > 1.1 Mq- The bright- 
est CMD-selected stars have V = 8.2 ( My ~ 0.2) or M ~ 2.5- 
3 Mq. The integration of the Kroupa (2001 ) IMF within these 



limits shows that we could miss nearly 90% of the cluster po- 
tential mass, even if, as discussed in Sect. 3, variations in the 
mass function just below our cut may significantly change the 
number of unstudied members. 

Despite this low percentage in mass of selected members, 
in Fig. g (left), two remarkable tails extend up to 1.5° from the 
cluster center, nearly perpendicularly to the Galactic Plane (the 
arrow points towards Z = 0). It is interesting to note the almost 
perfect circular geometry of the central parts of NGC 2516, 
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Fig. 3. Left: density (Log scale) of probable NGC 25 16 members, using star-counts on a B-R pair of Schmidt plates. The dashed 



circle represents the cluster diameter according to Lynga (1987). Dashed arrow is the proper motion (PM) of the cluster (cf. 
Table |3|) whereas the solid arrow points _L to the Galactic Plane towards Z = 0. Right: the colour-magnitude signal to noise 
(relative to the background, see text) with the mask (contour) used to select only high (bright regions) S/N areas. 





Fig. 4. Density of NGC 2287 stars (left) and pollution map for the same field (right). Definition of the arrows and dashed circle 
is the same as in Fig. |3|, with twice the NGC 2516 proper motion (PM) scale. Right: the IRAS 100 absorption is in greyscale. 
Small circles represent galaxies, as extracted from the SExtractor catalogue using a separation in a Log(area) vs. magnitude 
diagram. Many spurious extended source detections occur near the E border of the plate (X ~ —80' from the cluster center). 



despite the ideal location (£ ~ 270°) of this OC for observation 
of the predicted flattening by the radial galactic force gradient 
(see the discussion in Sect. 5). 

Except for a strong emission point north of the cluster (cor- 
responding to the "Toby jug" nebula IC 2220), the extinction 
due to dust is homogeneous, and no clumps in background 
galaxies correlate with the position of tails. Hence, no obser- 
vational bias can account for the presence of these huge tidal 
extensions. 



4.3. NGC 2548 (= M48) 

Despite being a bright cluster, NGC 2548 (later identified as 
the entry #48 in the Messier's catalogue) has not been studied 
recently. Here, we have used two pairs of blue and red survey 
plates, in order to check against possible variations in the sen- 
sitivity of the photographic emulsion which may bias the star- 
counts. No such effect was found. The OC is situated near the 
E border of the first pair (ESO/SRC 775) and on the W bor- 
der of the 2nd (ESO/SRC 776). Despite these locations on the 
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Fig. 5. Field centered on NGC 2548, two pairs of B & R plates. Left: the open cluster isodensity levels on the SRC775 plates. 
Right: the center of NGC 2548 lies only 20' East of the border of the SRC776 plate. Despite this, several features are nearly 
identical to the SRC775 plates (double core of the cluster, prolate shape and a tidal tail with a clump). See sect. 5 for the 
discussion about the clump of star A. Proper motion of the OC (PM), 4 times enlarged, is very uncertain. 



plates, the tidal tail extraction has proved consistent as shown 
on Fig. g. 

For the two pairs of plates including NGC 2548, we have 
chosen a strong cut, keeping only stars brighter than the in- 
strumental magnitude 13.0 (corresponding to about Z? cut ~ 14.8 
using a rough calibration based on GSPC-I standard stars), or 
Mb ~ 5.3. Again, we keep only stars more massive than about 
1-1.1 !Mq. Interestingly, the inner isodensity levels show a sec- 
ondary peak (double core). The projected axis ratio is around 
1.7:1 .0. The absorption in the region of the OC is very low 
(see Table ^) and uniform. Numerous extended sources were 
detected but no cluster of galaxies can account for the overall 
prolate shape and the SW tidal extension. 



5. Discussion 



It has been shown as well (Leon 1998) that the disk 



Theoretical considerations (Wielen 1975 ) have shown that the 
haloes of open clusters are flattened by the galactic gravi- 
tational field. Indeed, the gradient of the gravitational force 
in the galactic radial direction is stretching the OC: simula- 
tions (Terlevich 1987; Portegies Zwart et al. 2001) confirm 



that OCs take the shape of an ellipsoid with axes in a ratio 
X : Z ~ 2.0 : 1.0. The longest axis always points towards the 
Galactic Center whereas the smaller one is in the Z direction. 
This compressed shape has been observed only in a few clus- 
ters. In the Hyades, a flattening along the Galactic Plane was 
first suspected by Oort (1979 1. Perryman et al. (1998) recently 
confirmed, in their exhaustive study of the 3-D structure of this 
fundamental "calibrator" cluster, the prolate shape of the halo, 
which is well elongated towards the Bulge. In the Pleiades, a 
similar elongated shape roughly in the I direction was detected 
by van Leeuwen ( 1983 ) and definitively established by Raboud 
& Mermilliod (|l998a|). 



shocking compresses strongly the stellar cluster during the 
crossing of the galactic Plane. In the case of OCs, the heating is 
very efficient from the adiabatic component which includes the 
short period stellar orbits relative to the crossing time through 
the Plane. Figure gshows the contribution of the adiabatic heat- 
ing which is maximum for typical open clus ter ce ntral densities 
< 100 Mq pc~ 3 (e.g., Binney & Tremaine |l987| ). 

To quantify the effects of the vertical and radial (parallel to 
the Disk plane) tidal forces, we show in Fig. [7] the ratio of these 
two tidal components for a rather rich OC having a r t = 10 pc 
tidal radius. The galactic models are from Combes et al. (1999): 



the main difference is the scale height of the Disk which is set 
to 1 kpc for the Gal— 1 model and 400 pc for the Gal-2 model, 
leading to a major contribution of the vertical tidal force in the 
latter case owing to a larger gradient due to the thinner disk. 

However, the vertical component of the static tidal force is 
less important than the radial one, except on top of the Bulge 
where the radial component, in our model, is very weak. For 
a typical open cluster Z range, the vertical/radial ratio of the 
tidal force components is between 1% to 10%. Nevertheless 
the shock due to the time-dependent tidal force will induce a 
major heating of the OC crossing the Plane. Interestingly the 
vertical tidal force is stronger at very low-Z, where the open 
cluster population is mainly located, to decrease to a minimum, 
for a constant galactic radius, at Z ~ 150 pc (Gal-2 model) and 
increases towards higher Z where the vertical component from 
the Bulge and the whole Disk takes over. 

Hence, in addition to the static tidal field, all OCs suffer re- 
peated disk-shockings with the galactic Plane. Indeed one typ- 
ical open cluster at the solar radius crosses the Plane about 10- 
20 times before its dissolution in the Galaxy ( see be low), only 
leaving an OC remnant (de la Fuente Marcos 1998). For most 
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Fig. 7. Ratio (Log scale) of the Z-tidal force to the radial tidal force, parallel to the Galactic Plane for an open cluster having a 
tidal radius of r t = 10 pc. In a) we use the Gal-1 model (see Combes et al. 1999 for details) and b) the Gal-2 model. The main 
difference between the two models is the disk scale height, with 1 kpc for Gal-1 model and 400 pc for Gal-2 model. We did not 
take into account the dark matter halo. 
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Fig. 6. Contribution of the adiabatic heating during the disk 
cr ossing , usin g parti cle perturbation simulations (Leeuwin et 
al. 



1993 



Leon 



1998 ). The crosses represent the gain of energy 
?>E relative to the total energy Eo of the cluster. The dashed 
lines represent the dispersion in ?>E from different sets of sim- 
ulations. We indicate the range of rich OC central densities. 



of existing clusters, we expect to observe clumps of stars dis- 
rupted from the disk shocking, in the direction perpendicular to 
the Galactic Plane. These clumps should have a velocity (rela- 
tive to the OC center of mass) close to the velocity dispersion 
in the cluster. In an OC of ^500 SWq like the Hyades (evolved 



Table 4. Galactic positions and heliocentric space motions for 
the 3 selected open clusters, data coming from the compilation 
by Piatti et al. ( |1995 ). R p (resp. R a ) is the radius at the periph- 
elion (resp. at the aphelion) of the cluster orbit. 



Object 






Z 


# T 


U 


V 


W 


NGC 


kpc 


kpc 


pc 


pc 


km/s 


km/s 


km/s 


2287 


7.95 


8.97 


-130 


870 


+25 


+9 


-80 


2516 


5.60 


8.54 


-120 


200 


+2 


-24 


+5 


2548 


7.43 


8.99 


+ 160 


420 


+37 


+5 


-43 



T H= I 

Zmax | is the maximum height that the open cluster is sup- 
posed to reach above (or under) the Galactic Plane. 



intermediate-age cluster, like NGC 2548), the veloc ity disper- 
sion is about 0.2-0.4 km s~ 1 (de Bruijne et al. 200 1[ ). 

We can interpret the SW clump (noted A in Fig. around 
NGC 2548, in the direction of the Plane, as a remnant of the 
last disk-shocking. Given the distance of the clump to the cen- 
ter of NGC 2548 of about ^8 pc and taking a typical velocity 
dispersion of 0.3 kms -1 , we estimate the last shock finished 
ss26 Myr ago. On the other hand, following Dachs & Kabus 
(1989), the period Pz of the Z-oscillations of an open cluster is 
given by: 



Pz = 2k, 



(7) 



where Kz is the galactic vertical acceleration at the maximum 
height Z max the OC can reach. We use the acceleration values 



given by Vergely et al. ( 2001 ), based on recent estimates of the 
disk scale height D = 240 pc and surface mass density Eq = 
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48 5Wq pc 2 , whereas the local mass density is taken as po = 
0.076 pc" 3 . 

The results of the compilation of the absolute (U, V, W) mo- 
tions of the three clusters by Piatti et al. (1995) are shown in 
Table [|. In the case of NCG 2548 however, different values 
were found recently by Wu et al. (2001). These authors de- 
duced (U,V, W) = (l,221,3)kms from an exhaustive proper 
motion survey, with H = |Z max | = 170 pc (taking the nearby 
distance estimate of Jm48 — 530 pc, see Clana 1985). 



Using the W and H values of Piatti et al. (|1995J), we can es- 
timate that the last crossing of NGC 2548 through the Galactic 
Plane finis hed ab out ^40 Myr ago, whereas taking the data of 
Wu et al. ( |200l[ ) the last shock occurred ~20 Myr ago. Both 
results are consistent with our dating derived from the clump. 
We note that with an age of about 360 Myr (see Table |), NGC 
2548 is a dynamically evolved OC which has already suffered 
about 7 or 8 disk-shockings. 

Combes et al 



It has also been shown (Leon 1998 



1999) 



that during the disk-shocking the shape of the cluster will 
change rapidly before the relaxation, which is never reached 
for an OC constantly oscillating in the Plane. The difference in 
shape between NGC 2516 (round core) and both NGC 2287 
& NGC 2548 (pronounced prolate shape) can be accounted 
for a transient phase during the shock. It could also be due 
to effects of encounters with passing molecular clouds which 
can reduce, or randomize the orientation of the flattening (see 
Theuns|l99lal[l992bi. 



In order to estimate the importance of the disk-shocking 
in the evolution of the OCs, we have performed A^-body sim- 
ulations. We us ed the so-called perturbation particles method 
(Leeuwin et al. 1993| ). To solve the collisionless Boltzmann 



equation, this technique combines an analytical description of 
the equilibrium state and a numerical evaluation of the lo- 
cal density perturbation. Indeed perturbation particles are well 
suited for a collisionless system, and can be used for simulat- 
ing a disk-shocking because the relaxation time of the cluster is 
still larger than the crossing time. The simulated cluster is mod- 
eled by a Plummer's sphere for which the gravitational poten- 
tial a nd the distribution function are known analytically (e.g., 
Leon|l998|). 



We stress that our simulations are only aimed at giving an 
estimation of the mass loss. In particular, we have not included 
the effect of binarity and mass segregation. Binary stars, pri- 
mordial or exchanged, seem to be less important for the evolu- 
tion of rich clusters like those considered in our models (see de 



la Fuente Marcos 1996) and may be regarded as more massive 
stars, hence defering the problem to that of mass segregation. 

Furthermore, the perturbation particles method does not 
take into account the mass segregation process, but during one 
crossing this latter will be slow (see Combes et al. 1999). The 
spatial mass function is certainly different between the core and 
the outskirts of the cluster (e.g., Raboud & Mermilliod 1998b). 
However, during a crossing the tidal effects are independent of 
the stellar masses. As a result, ignoring the mass segregation 
is not expected to strongly bias the total mass of the tidal tail 
(only the mass function will be different, with an excess of low 
mass stars). 




Fig. 8. Z position of the open cluster vs. the time (expressed in 
Af-body time units) in the perturbation particle simulation. 




Fig. 9. Mass loss relative to the total mass of the cluster vs. 
time (A^-body units) in the perturbation particle simulation. 



We have followed an OC starting at 600 pc (see Fig. ||) 
above the Plane at the solar galactic radius with a total mass of 
1000 5Wq and a core radius r c of 1 pc. We have chosen such 
a high-Z cluster to avoid strong gradient at the begi nning of 



the simulation. Old "survivors" like M67 (Fan et al. 1996) or 



NGC 188 (Sarajedini et al. |1999| ) have total masses of t his or - 
der, and follow such high-Z orbits (see Carrara & Chiosi 1994 ). 
Figure ^| shows that the OC mass loss by the disk-shocking is 
about 1.5% of its total mass in one Z-oscillation. It can be no- 
ticed that some stars are re-bound to the cluster after the cross- 
ing. 

Despite the lower intensity of the static Z-tidal force rela- 
tive to the radial one, the time -dependent Z-tidal force during 
the crossing strongly speeds up the destruction of the OC. The 
global contribution of the disk-shocking can be estimated to be 
about 10 to 20% for such a high-Z cluster, and likely more for a 
younger open cluster which experiences stronger Z-force gra- 
dients in the vicinity of the Plane (see above). We can try to 
relate this mass-loss to the mean life duration of open clusters 
in the Galaxy. 
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Using a pre-release of the catalogue of Loktin et al. (2001) 
available on the Base des Amas web database, we show in 
Fig. |To| the age distribution of the 423 open clusters in their 
sample. As mentioned above, OCs suffer different destructive 
processes (disk-shockings, encounters with GMCs, tidal field 
harassment) which lead to the rapid depletion of the whole pop- 
ulation. It appears that the population age, excluding both very 
young (< 10 yr) and very old (> 1 Gyr) OCs, can be mod- 
eled by an exponential law of the following form (see the fit in 
Fig.©: 



N = N*e~ 



(8) 



where T+ is a time-scale for OCs (primarily in the solar 
neighbourhood, as catalogues become incomplete for distant 
clusters) which corresponds to their destruction time-scale. 
Assuming a constant star-formation rate over the last 10 9 years 
in the solar neighbourhood (e.g., Rocha-Pinto et al. 2000), 
the OC destruction time-scale is of the order of 6 x 10 s yr, 
somewhat superior to previous observational estimates (Lynga 



1982; Janes & Phelps 1994). Selection effects in different 



Galactic OCs catalogues may explain this small discrepancy. 
Interestingly, the value of 600 Myr is in good agreement with 
the OC disintegration time in the field x , defined by Rroupa 
( 1995 ) in its realistic simulations. 



As stated before, in the solar neighbourhood, disk shocking 
seems to be very important in the destruction of OCs. An OC 
at the solar radius suffers about ten to twenty crossings dur- 
ing this time-scale, which leads to an upper-limit for the mass 
loss efficiency during one disk-crossing between 5 to 10%. This 
value is clearly more important than for globular clusters (Leon 



1998) which are much more massive and concentrated stellar 



aggregates. 

6. Conclusions 

We have detected important tidal extensions around 3 bright 
open clusters using photometric selection on wide-field pho- 
tographic plates. Density isophotes of cluster probable mem- 
bers go further away from the ordinary limits given for these 
objects. In two cases (NGC 2287 & NGC 2548) we have ob- 
served strong flattening of the cluster that we interpret as a 
consequence of the disk-shocking. We stress that the adiabatic 
heating is important for the open cluster evolution, as well as 
the time varying galactic field. In spite of the weakness of the 
static Z-tidal force field, the time varying tidal field during the 
disk-shocking will indeed affect deeply the dynamical evolu- 
tion of OCs. Its effects, computed with the perturbation parti- 
cles method, account for about 20% of the mass loss. We also 
estimate the destruction time-scale of open clusters in the so- 
lar neighbourhood to 600 Myr. In the case of NGC 2548, it 
has been possible to estimate that the last shocking with the 
Galactic Plane occurred about 30 Myr ago. This value should 
be confirmed by a complete proper motion and radial velocities 
study. Advent of super-MOS (e.g. Flames) and, at medium- 
term, Gaia will provide hundreds of V r and very accurate 
proper motions even for low-mass members of known OCs. 
Moreover, accurate multi-band, deep wide-field CCD photom- 
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Fig. 10. Distribution of the numbe r of O Cs vs. the age (Log) 
from the catalogue of Loktin et al. ( 2001 ). The exponential fit, 
excluding the youngest and the oldest Galactic open clusters, 
with a time-scale of 612 Myr is overplotted as a dashed line. 

etry of a sample of open clusters will provide much better selec- 
tion of members, including sub-solar mass stars. The compar- 
ison of these new observations with more realistic simulations 
will be the subject of forthcoming papers. 
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